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Abstract. Recent observations provide determinations of 
individual masses, chemical composition and metallicity of 
the components of the spectroscopic and interferometric 
binary iPeg (Boden et al. 1998). Using updated physics, 



al. 1994), stellar diameters (e.g. Koechlin & Rabbia 1985 
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to calibrate the system, we have computed using the stel- 



lar evolutionary code CESAM (Morel 1997), evolution 



ary sequences of stellar models with the masses of iPegA 
1.326 M Q and tPegB 0.819 Mq (Boden et al. loc.cit) and 
with different values of the mixing- length parameter a, the 
helium Y and the heavy element Z initial mass fraction 
with the constraint of the observed metallicity. Adopting 
effective temperatures and luminosities, as derived from 
observations with the bolometric corrections, and the em- 



pirical scale of temperatures of Alonso et al. (199J, |l996j) , 
we find a A = 1.46, a B = 1.36, Y = 0.278, Z = 0.017. 
The evolution time, including pre-main sequence, is found 
within ~ 40 My <; t cv < 0.5 Gy. The calibrated models 
of tPeg.A and B are non homogeneous zero age main 
sequence stars with the evolutionary time t QV = 56 My. 
Due to the large uncertainties of their determinations, 
the values derived for the mixing-length parameters are 
smaller than the solar one but however marginally com- 
patible with it. Our results ought to be improved as soon 
as a more accurate value of the magnitude difference in the 
V filter will be available. Detailed spectroscopic analysis 
for both components looks practicable, so it is urgently 
needed. 

Key words: Stars: binaries: spectroscopic - Stars: bina- 
ries: visual - Stars: fundamental parameters - Stars: indi- 
vidual: l Peg 



1. Introduction 

Along the last decade, owing to significant improve- 
ments of our knowledge of fundamental stellar data, 
trigonometric parallaxes with HIPPARCOS (ESA |1997| ) , 
spectroscopic-speckle mass determination (e.g. Scarfe et 



DiBcnctto & Bonncau 1991), photometric (Mermillod et 
al. 1997) and spectroscopic data (Cayrel de Strobel et al. 
1997), the constraints on stellar modeling became more 
stringent leading to improvements of stellar evolution the- 
ory. Meanwhile, physics, e.g. opacities (Iglesias & Rogers 
1996), equation of state (Rogers et al. [1996 ; Dappen 1996), 
thermonuclear reaction rates (Adelberger et al. 1995 , An- 
gulo et al. 1999), have been improved. Most of these new 



data are nowadays available in free access in data banks. 

With the exception of the Sun, modeling a single star 
is not a closed problem because the number of indeter- 
minate parameters is larger than the observed ones. For 
some well known binar ies, e.g. a Cen (Demarque et al. 



1986; Noels et al. 1991), with the reasonable hypothesis 
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of a common origin for both components, i.e. same initial 
chemical composition and age, it is possible to confront ob- 
servations and theoretical stellar evolution sequences. One 
can derive estimates for age, helium content and metallic- 
ity, fundamental quantities for our understanding of the 
galactic chemical evolution. One derives also the mixing- 
length parameter a - also named "convection parameter" . 

For a fixed physics, the value of the mixing-length pa- 
rameter 1.0 ^ q q 2.5 derived from a calibrated solar 
model is currently used to model other stars. But, one can 
ask why the value adjusted for the present Sun should 
be applied to stars with different masses, ages and/or 
chemical composition. Fixing the physics, Fernandes et 
al. (1998) have shown that a Q is valuable for a sample of 
low mass binaries 0.6 M© < M* 1.0M Q in the metal- 
licity range -0.65 < [F c/H] < +0.0. Indeed, for a Cen 
Fernandes & Neuforge ( |1995| ) obtained a value close to 
the solar one. More recently, Pourbaix et al. (1999) using 
improved mass and gravity determinations for both com- 
ponents of a Cen were able to relax the assumption of a 
unique value for the convection parameter. They obtained 
a A = 1.86±0.80 (resp. a B = 2.1±0.80) for the A (resp. B) 
component. Due to the range of uncertainties, this result 
does not rule out the assumption of a constant value for 
the mixing-length parameter for masses and metallicities 
smaller or slightly larger to solar values. 
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Up to now, the calibrations of double stars have espe- 
cially concerned main sequence stars with ages of few Gy. 
As soon as new data become available, it is important to 
extend the sample of calibrated binaries. In most cases 
the individual masses are inaccurately known and the 
isochrones are not really constrained. The spectroscopic 
and interferometric binary t Peg. is a particularly inter- 
esting object since the two masses are accurately known 
(Boden et al. loc. tit.). Moreover, the mass and the metal- 
licity of the primary are slightly larger than solar. It is a 
"solar like" star whose seis molog ical prop erties have been 
stu died (e.g., Houdek et al. |l995| , Houdek |l996| Michel et 



al. |l999). Some stars of this group are COROT targets 



(Baglin 1998 ). The values of the mixing- length parame- 
ters, as derived from calibration for iPeg.A and B, will 
provide a test for the 2-D numerical calculations of com- 
pressible convection (Ludwig et al. 1999 ). 

Some authors have already derived estimates for the 
age of the iPeg. binary system. A value of <~ 80 My 
was derive d fro m the obse rved lithium depict ion ( Fekel 
& Tomkin |l983| , Gray |l984|) . Lyubimkov et al. |l99l| have 
obtained the slightly larger value of 170 ± 80 My, from 
comparison with evolutionary models. None of these es- 
timates is based on evolutionary models taking explicitly 
into account, as we do, the microscopic diffusion and grav- 
itational settling, basic phenomena in stellar modeling. 

In this paper we report an attempt to calibrate tPeg. 
using stellar evolution models including pre-main sequence 
and microscopic diffusion. It is organized as follows: In 
Sect. | we collect and discuss the observations. Section || 
is devoted to derivations of effective temperatures and 
luminosities. The method of calibration is described in 
Sect. |]. Section || presents the stellar modeling procedure. 
In Sect. § we give the results with emphasis on a seismo- 
logical analysis. We discuss the results and conclude in 
Sect. 

2. Observations of i Peg. 

iPegasi (HR 8430, HD 210027) is a nearby (d ~ 11.75 pc), 
a wonderful double lined spectroscopic and interferometric 
double star. The binar ity wa s discovered more than a cen- 
tury ago by Campbell (1899), who notices "The velocity of 
6 Pegasi in the line of sight is variable" . Table [l] shows the 
astrometric, photometric and spectroscopic observational 
data retained for the present investigation. 



2.1. Orbit, masses, parallax and diameters 

Recently, Boden et al. {loc. cit.) using observations of the 
Palomar Testbed Interferometer computed the visual or- 
bit of l Peg. and derived the individual masses. The orbital 
period is accurately known p = 10.213033 ± 1.3 10~ 5 day 



Table 1. Astrometric, photometric and spectroscopic 
adopted data of i Peg. my, Amy (resp. tok, Attzk) are the 
global apparent magnitudes and magnitude differences in 
the V (resp. K) band. w OI h is the orbital parallax. Ma, 
Mb are the masses, [Li]a, [Li]b the lithium abundances 
(H=12), of A and B components, respectively. [Fe/H] is 
the metallicity. The flag signals a value estimated by 
ourselves. 



Spectra 

mv 

Amv 

m K 

AmK 

■Ci7 or b 

Ma 
Mb 

Ma 
Mb 
[Fe/H] 



F5V - G8V 
3.771 ± (0.001)* 
2.68 ± (0.05)* 
2.656 ±0.002 
1.610 ±0.021 
86.91 ± 1.0 mas 
1.326 ±0.016 M 
0.819 ± 0.009 M Q 
2.9 ± (0.4)* dex 
2.6 ± 0.6 dex 
(±0.00 ±0.1)* dex 



Fekel & Tomkin ( 1983 ) 
Fekel & Tomkin ( fl983| ) 
See text 



Fekel & Tomkin ( 1983 ) 
Mermillod et al. (1997) 
Fekel & Tomkin (p83p 
Bouchet et al . (|199 
Boden et al. (|1998|) 
Boden et al. ( 1998 ) 
Boden et al. ( 1998 ) 
Boden et al. (1998) 



(Fekel & Tomkin 1983). The spectroscopic and visual or- 
bits for tPeg. are each statically consistent with being 
circular. That corroborates the result of Zahn & Bouchet 



( |1989D who predict, for close binary systems with masses 
ranging from 0.5 to 1.25 M@, that the orbit is circularized 
during the Hayashi phase. Moreover, since the period of 
l Peg. is larger than w 8 days, it is expected that the main 
sequence is reached with components rotating faster than 
the orbital rate, Zahn & Bouchet (loc. cit.). 

The orbital parallax derived by Boden et al. (loc. cit.), 
as seen in Table [j], is in close agreement with the trigono- 
metric tutrig = 85.06±0.71mas (ESA [l997 ) determination 
by HIPPARCOS. Though the accuracy of the HIPPAR- 
COS parallax is slightly higher, we have used the Boden 
et al. value for global consistency. Boden et al. (loc. cit.), 
assuming two different limb darkening laws, have obtained 
consistent estimates of diameters for each component: 
© A = 0.99 ±0.05 mas, © B = 0.70 ±0.10 mas. The illumi- 
nating Fig. 1 of Boden et al. (loc. cit.) shows that the pre- 
vious values of diameters are not large enough for eclips- 
ing, despite the high inclination i — 95.85° ±0.22° and the 
small semi-major axis a — 10.33±0.10 mas. From the pre- 
vious interferometric data and orbital parallax one derives 
the "interferometric" radii of i? A in = 1-224 ± 0.025 R & 
and i?Bin = 0.872 ± 0.050 Rq, respectively, for primary 
and secondary. As emphasized by Boden et al. (loc. cit.), 
the sum of the radii is consistent with the absence of evi- 
dence of photometric eclipse. For the gravities one obtains: 
log 10 3Ai „ = 4.39 ± 0.23 and log 10 g Bin = 4.47 ± 0.65. 



2.2. Spectroscopy 

From analyses of 2.5 A per mm CCD spectra, Lyubimkov 



et al. ( |199l|) have derived effective temperatures T e g Asp = 
6750 ± 150 K, T cffBs p = 5350 ± 350 K, and gravities 
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logic 9 As 



4.35 ± 0.05, log 10 5 bs 



4.57 ± 0.10. 



With the values of the masses derived by Boden et al. 
(Zoc. cit.), that correspond to a "spectroscopic" luminosity 



of L Asp = 3.024±0.62£ Q and L Bsp = 0.444 ± 0.219 L , 
respectively for the primary and secondary, and to radii 
i? Asp = 1-274 ± 0.07 R e and R B 



sp 



0.772 ± 0.089 R 



©• 



Another effective temperature value T e g = 6488 K, de- 
rived from spectroscopic analysis, was just recently listed 
by Boesgaard et al. (199£). Note that the constraint on 



the diameters, due to the absence of evidence of photo- 
metric eclipse (Boden et al. loccit.), is also fulfilled with 
the "spectroscopic" values of the radii. 

2.3. Chemical composition 



Lyubimkov et al. (1991) have qualified the chemical com- 
position as "close to normal". The Catalog of [Fc/H] 
(Cayrel de Strobel |l997| ) recommends [§] = +0.10. An- 
other value of the metallicity [^] = —0.13 is due to 



Duncan ( |l98l| ). According to Boesgaard et al. (1999) the 
metallicity is [^] = —0.08. Faced with these rather scat- 
tered data, we have adopted for our calculations = 
0.00 ± 0.10 (Cayrel de Strobel, private communication). 
Except for lithium (see next), none of the authors has 
reported significant differences between the surface chem- 
ical composition of primary and secondary. It is an indi- 
cation that both components of tPeg. are newly formed 
stars, since the gravitational settling, which acts more effi- 
ciently in a 1.326 M than in a 0.819 M , still has not have 
enough time to differentiate significantly the surface chem- 
ical composition of primary and secondary (e.g. Morel & 
Baglin [l999[ ). Therefore, the initial metallicity of tPeg. is 
certainly very close to the nowadays observed value. 

Table |l| shows the adopted values of lithium abun- 
dances. They are in agreement with the measurements 
of Lyubimkov et al. ( 1991 ): [Li] a = 3.25 for primary and 
[Li]b =2.58 for secondary. According to the standard evo- 
lution theory, for an isolated star of solar metallicity and 
less massive than M* Ss 0.85 M , all lithium is destroyed 
at the end of the fully convective phase of the pre-main 
sequence. That results from how fast the convective zone 
recedes from the core where the lithium is depleted - recall 
that lithium burns as soon as the temperature is greater 
than T 2; 2.4 MK. In the core of a late type star, the in- 
crease of temperature is slow, so are the decreases of opac- 
ity and radiative gradient. Then, the fully convective state 
lasts longer than in a more massive star. With respect to 
an isolated star, tPeg. B is lithium under-depleted. For 
not isolated stars under-depletion can occur on the main- 
sequence in binary systems which are tidally locked for 
orbital periods below ~ 8 days (Zahn 1994). The period 
of tPeg. (10 days) is slightly larger than this limit, as 
supported by the no-synchronous rotation found by Fekel 
& Tomkin ( |1983| ). But, Gray ( |1984[ ) deduced, from the 
analysis of higher resolution spectra, that primary and 
secondary are slow rotators tAsini = 6.5 ± 0.3kms _1 , 



vb sin i = 5. ± 1 kms^ 1 , and argued that both components 
are in synchronous rotation. As previously conjectured, 
the components may have reached the main sequence with 
rotation velocities larger than the orbital one and the indi- 
vidual rotation velocities may have decreased, the system 
being again tidally locked. Figure [l] shows that the locus 
of l Peg. B in the HR diagram excludes the possibility that 
the secondary is nowadays still in the pre-main sequence, 
before the lithium burning phase which occurs soon after 
the end of Hayashi's track. Therefore, the under-depletion 
of lithium, observed in 6 Peg. B, is unexplained. There is 
an alternative either, i) along the pre-main sequence, for 
some unknown reason, the mixing of the convection zone 
is slowed and lithium is not destroyed or, ii) at the be- 
ginning of the main sequence, coming from somewhere, 
lithium has been provided to the outer convection zone 
of the secondary. Due to the distance between the compo- 
nents a w 0.12 AU there is no possibility for mass exchange 
between the components (F. van't Veer, private commu- 
nication). Then the first part (i) of the alternative appears 
to be the most reliable. That reinforces the conjecture that 
l Peg. is certainly composed of two newly formed stars. In 
our calculations, to mimic the unknown process leading to 
the under-depletion of lithium in the secondary, the con- 
vective zones will be only mildly mixed in the course of 



the pre-main sequence (see Sect. 5.2). 



2.4- Photometry 

i Peg. A is a MK standard and its photometry is accu- 
rately known. Another value of the K apparent magnitude 
m K = 2.623 ± 0.016 is due to Carrasco et al. (|199l[). For 



the magnitude difference Amy in the V band we have 
adopted the estimate of Fekel & Tomkin ( 1983| ). Their 
analysis is based on spectroscopic and photometric argu- 
ments. Comparing the strength of iron lines of t Peg. B 
and e Eri, a K2V star, they derive, for t Peg. an estimate 
of the magnitude difference ArriR = 2.5 ± 0.3 in the R 
band. Such a value of Ampt characterizes stars between 
G5V and K0V. Then, they deduce Amy = 2.68 from the 
sensitivity of the global K magnitude for binary systems 
composed of a F5V primary and secondaries with spectral 
types between G5V and K0V. Fekel & Tomkin (loc. cit.) 
do not give the uncertainty 5 Amy of their estimation of 
Amy; it is difficult to derive such value from their subtle 
analysis. We have estimated SAmy = ±0.05. This quan- 
tity is the most poorly known parameter of our calibra- 
tion. Lyubimkov et al. ( 1991 ) have reported the estimate 
Amy = 2.13 which significantly differs from the value de- 
rived by Fekel & Tomkin (loc. cit.). 

3. Luminosities and effective temperatures 

Faced with the scattered values of luminosities and ef- 
fective temperatures, derived eith er f rom spectroscopic 
(Sect. |2.2[) or photometric (Sect. 2A) observations, we 
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Table 2. Physical parameters derived from data of Ta- 
ble using n7 or b- Bck is the bolomctric correction relative 
to the K magnitude. Mboi, My (resp. Mr) are respec- 
tively the bolometric and absolute magnitudes in the V 
(resp. K) band. 



Parameter 


l Peg. A 


u Peg. B 


Bck 


+0.97 


+1.82 


Mboi 


3.54 ± 0.08 


6.01 + 0.09 


M v 


3.56 ± 0.03 


6.24 + 0.07 


M K 


2.57 ±0.03 


4.18 + 0.03 


Teff 


6642 ± 52 K 


4992 ± 115 K 


L/Lq 


3.041 ± 0.23 


0.3145 ± 0.027 



make the choice to use the photometric data alone. Ta- 
ble |^ gives the physical parameters of primary and sec- 
ondary as derived from data of Table [l]. The bolometric 
corrections for the K band are obtained using the bolo- 
metric corrections of Alonso et al. ( [1995P with a correc- 
tion of the zero point giving a bolometric magnitude of 
4.75 to the Sun (Cayrel, private communication). The un- 
certainty of the method is of the order of w ±0.05 mag- 
nitude. The e ffectiv e temperatures are derived from the 
Alonso et al. ( 1996 ) empirical scale of temperature, valid 
for the low main sequence from F0V to K0V. The uncer- 
tainty of the fit is of the order of ST c g ps ±30 K. We are 
aware that many tidally locked binary systems, having a 
lithium under-depleted late type component, are recog- 
nized as chromospherically active (Barrado y Navascues 
1997). As iPeg. is not identified as chromospheri- 



et al. 

cally active (Barrad o y Navasc ues et al. loc. cit.) the use 
of the Alonso et al. ( 1995 , 1996) procedures is valid. From 
the previous photometric and orbital data one derives the 
"photometric" radii of i?Aph = 1.319 ± 0.023i?o for pri- 
mary and i?e ph = 0.751 ± O.O27i?0 for secondary. Indeed, 
for the gravities one obtains: log 10 3a ph = 4.32 ± 0.02, 
log 10 3b ph = 4.60 ± 0.03, values in agreement with their 
"spectroscopic" determination of Sect. 2.2. Note that the 
constraint on the diameters, due to the absence of evi- 
dence of photometric eclipse (Boden et al. loc cit.), is also 
fulfilled with the "photometric" values of radii. 

The effective temperatures, luminosities and radii of 
F5V and G8V standard stars are respectively 6440 K, 
3.20£ e , 1.44 R Q and 5570 K, 0.66 L , 0.874 R Q (Schmidt- 



Kaler 1982). According to radii and luminosities, derived 



previously either from photometric data or from spectro- 



scopic data (Sect. 2.2), both components of tPeg. are un- 
der luminous and have smaller diameters than standards. 
That reinforces the conjecture of a small age. tPeg. being 



certainly a young binary system, the models have to take 
properly into account the pre-main sequence evolution. 

Uncertainty domains. For each component the luminos- 
ity and the effective temperature, as derived from Alonso 
et al. (loc. cit.) semi empirical adjustments, are functions 
of absolute magnitudes My and Mk . Then, the luminos- 
ity and the effective temperature have the following non- 
linear functional dependency with respect to the observed 
parameters: 

L = L(my, Amy, tuk, AmK, w, [Fe/H]), 
T cf f = T of f (m v , Amy, m K ,Am K ,n7, [Fe/H]). 

Only the uncertainties on the observed magnitudes, the 
distance and the color index can be assumed to be Gaus- 
sian. The statistical behavior of the metallicity deter- 
mination, of the estimated V magnitude difference (see 
Sect. 2.4) and of the fits of Alonso et al. (loc. cit.) are 
unknown. Therefore, due to the non-linear dependence of 
the effective temperature with respect to variables some 
of which are correlated and not Gaussian, we have esti- 
mated the uncertainty on the effective temperature using 
the sum: 





dT oB 


Smy H 








dmy 




d Amy 


+ 


dT cS 
<9tok 


Sm K + 


dT cS 
<9Atok 


+ 


dT cS 


5za + 


dT eS 


dm 


d [Fe/H] 



AT c g « - c Smy + - - c — 5 Amy + 

i5Atok + 
8[Fe/R] + ST^, 



instead of the standard square root of the sum of the 
square of the standard deviations. A similar relation holds 
for the luminosity. Here, 5x is the uncertainty on x. Table 
shows the uncertainties derived in such a way for effective 
temperatures and luminosities. They are maybe slightly 
optimistic but give a reasonable order of magnitude. 

4. The method of calibration 

Fixing the initial mass and the physics, the calculation of a 
stellar model needs 4 unknowns namely, the mixing-length 
parameter a, the helium To and heavy elements Zq initial 
mass fractions and the age. To model a binary, assuming a 
common origin to each component, the age and the initial 
chemical composition are identical for primary and sec- 
ondary. There are only 4x2 — 3 = 5 unknowns. There are 6 
observables: luminosity, effective temperature and metal- 
licity for each component. In case of i Peg., many presump- 
tions of a small age allow to assume that the observed 
surface metallicity (y) t p cg . is the same for both compo- 



nents and then close to its initial value (see Sect. 2.3). 
Therefore, there remain 5 observables and 5 unknowns: 
the mixing- length parameters a t p e g.A and a t p Gg . b of pri- 
mary and secondary, the initial helium Y and metal Z 
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contents and the age i t p cg .. As 1 = X + Y + Z, the initial 
values of helium and metals are related by: 

(i-y ) z 

This relationship holds only because the initial and present 
day mctallicities are alike for both components; otherwise 
Yq and Zq need to be derived independently. Therefore, 4 
calibration parameters remain: 

C*i,Peg.A, O^Peg.B, ^0) ^iPeg.- 

Since we have adopted ] t p og . = 0.0, we used (^■) i ,p cg . = 



(§) = 0.0245 (Grevesse & Noels. |1993| ) 



Search of a solution. In a first time we have explored sys- 
tematically the domain of the calibration parameters. In 
a second step, more or less empirically, we have refined 
the convection parameters and the initial helium mass 
fraction. Doing that we are guided by two facts, i) in 
the HR diagram an increase of the convection parame- 
ter moves the evolutionary track towards larger tempera- 
tures at about constant luminosity and, ii) an increase of 
the helium content shifts the evolutionary track towards 
larger effective temperatures and luminosities. For stars 
with convective envelope, the first point (i) results from 
the fact that an increase of the efficiency of the convection 
deepens the base of the convection zone and does not affect 
the core. The luminosity is then almost unchanged, while 
the adiabat of the external convection zone is adjusted 
on a larger temperature value at the base of the convec- 
tion zone, with the result of an increase of the effective 
temperature. The second point (ii) can be understood in 
the following way: an increase of the helium mass fraction 
leads to an increase of the mean molecular weight which, 
by turn, diminishes the pressure gradient. The quasi-static 
equilibrium is conserved via a temperature increase. A 
concomitant luminosity increase results from the enlarge- 
ment of the nuclear energy generation magnified by the 
large power law dependence of nuclear reaction rates, with 
respect to temperature. In summary, in the HR diagram 
a change of a translates the locus of a model in effective 
temperature at constant luminosity, while one can adjust 
changes of a and Y in order to translate the locus of a 
model in luminosity at constant effective temperature. 

Finally, the steepest descent algorithm (e.g. Conte & 
deBoor 1987, Noels et al. 1991) is used to improve the cal- 



ibration parameters. As this gradient method converges in 
the vicinity of stable solutions we have only consider as a 
"solution" any set of calibration parameters for which the 
gradient method converges. We have kept for the calibrated 
models the solution which best fits the central values of 
observational boxes in the HR diagram for each compo- 
nent. Then we have estimated the domains of variations 
allowed to the parameters in order that the loeii of the 
models remain in these observational boxes. 



5. The physics of the models 

Basically the p hysics used in the models is the same as 
in Morel et al. ( 1997 ). The standard assumptions of stel- 
lar modeling are made, i.e. spherical symmetry, no ro- 
tation, no magnetic field, no mass loss. Each evolution 
is initialized with a homogeneous zero-age pre-main se- 
quence model in quasi-static contraction with a tempera- 
ture at center close to the onset of the deuteron burning, 
i.e. T c ~ 0.5 MK. We shall call the "time of evolution" 
of a model the time t cv elapsed from initialization. In a 
sequence of models we shall designate by "model of zero 
age main-sequence" (ZAMS) the first model, if any, where 
the nuclear reactions dominate gravitation by more than 
99% as the primary energy source. Typically ZAMS oc- 
curs after ~ 40 My of evolution for t Peg. A and ~ 70 My 
for i Peg. B. 



5.1. Nuclear and diffusion network. 

The general nuclear network we used contains the fol- 
lowing species: 1 H, 2 H, 3 He, 4 He, 7 Li, 7 Be, 12 C, 13 C, 
14 N, 15 N, 16 0, 17 and Ex; Ex is an "Extra" fictitious 
mean non-CNO heavy element with atomic mass 28 and 
charge 13 (Ex ~ 28 Al) which complements the mixture 

i.e., Xex = 1 — zJ-=ih ^ i with Xi as the mass fraction 
of the species labeled with i = 1 H, . . . , 17 0. With respect 
to time, due to microscopic diffusion processes, the abun- 
dances of heavy elements are enhanced towards the cen- 
ter; Ex mimics that enhancement for the non CNO metals 
which contribute to changes of Z , then to opacity varia- 
tions but not to the nuclear burning and nuclear energy 
generation. We have taken into account the most impor- 
tant nuclear reactions of PP+CNO cycles (Clayton, 1968 ). 
The relevant nuclear reaction rates are taken from the 
NACRE compilation (Angulo et al. 1999| ) with the reac- 
tion 7 Be(e~, ^ 7) 7 Li taken from the compilation of Ad el- 
berger et al. ( 1998 ). Weak screening (Salpeter 1954 ) is 
assumed. Owing to the fact that the nowadays observed 
surface chemical composition is solar and close to the ini- 
tial mixture (see Sect. ||), the initial fractions between the 
heavy elements within Z are set to their photospheric so- 
lar present day values. More details about this procedure 
can be found in Morel & Baglin ( 1999[). We have used the 
meteoritic value (Grevesse & Sauval |1998| ) for the initial 
lithium abundance, [Li] = 3.31 ± 0.04 (H = 12). For the 
calculations of the depletion, the lithium is assumed to be 
in its most abundant isotope 7 Li form. The initial abun- 
dance of each isotope is derived from isotopic fractions and 
initial values of Y = 3 He + 4 He and Z in order to fulfill 
the basic relationship X + Y + Z = 1 with X = 1 R + 2 H. 
Microscopic diffusion is des cribed by the simplified formal- 
ism of Michaud & Proffitt ( 1993 ) with each heavy element 
as a trace element. 
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5. 2. Equation of state, opacities, convection and 
atmosphere. 

We have used the OPAL equation of state (Rogers et al 



1996 ) and the opacities of Iglesias & Ro gers (1996 ) for 
the solar mixture of Grevesse & Noels ( 1993| ) comple- 
mented, at low temperatures by Alexander & Ferguson 
(|1994D opacities. 

In the convection zones the temperature gradient is 
computed according to the standard mixing-length the- 
ory (Bohm-Vitensc 1958] ) . The mixing- length is defined as 



/ = aH p , where H p is the pressure scale height. The con- 
vection zones are mixed via a strong full mixing turbulent 
diffusion coefficient d{ m — 10 13 cm 2 s _1 which produces 
a homogeneous composition (Morel 1997). As quoted in 
Sect. ||, during the pre-main sequence, we used a soft mix- 
ing of the convection zones. We found empirically that 
the mixing generated by a turbulent diffusion coefficient 
drnm = 10 6 cm 2 s -1 provides lithium depletions close to 
the observed values for both components. We are aware 
that this procedure is not fully satisfactory: first, as al- 
ready quoted, we do not know why and how the lithium 
destruction is annihilated in l Peg. B during the pre-main 
sequence, second, we suppose a priori that the observed 
depletion corresponds to the onset of the main sequence 
and, third, we do not know how the unknown physical 
process which acts on lithium performs on other chemi- 
cals. 



Following the prescriptions of Schaller et al. ( 1992 ) the 
models take into account overshooting of mixed convective 
cores over the distance O v — 0.2 min(_ff p , R co ); R co is the 
radius of the convective core. 

A n atm osphere is restored using Hopf's T(r) law (Mi- 
halas 1978 ). The connection with the envelope is made at 
the Rosseland optical depth rb = 10, where the diffusion 
approximation for radiative transfer becomes valid (Morel 
et al. 1994). In the convective part of the atmosphere, a 
numerical trick (Henyey et al. 1965| ) is employed in connec- 
tion with the purely radiative Hopf law in order to ensure 
the continuity of gradients at the limit between the at- 
mosphere and the envelope. At each time step, the radius 
R+ of the model is taken at the optical depth ~ 0.6454 
where T(t±) = T g; the mass of the star M*, is the mass 
inside the sphere of radius i?*. The external boundary is 
located at the optical depth r ext = 10~ 4 , where the density 
is fixed to p(r C xt) = 3.55 10 _9 gcm~ 3 . 



Table 3. Characteristics of a calibration of the 6 Peg. sys- 
tem. i? cz and R co are respectively the radius of the base 
of the external convection zone and of convective core (in- 
cluding overshoot). At the center, T c , p c , X c , Y c are re- 
spectively the temperature (in MK), density (in gem -3 ), 
hydrogen and helium mass fractions. 





l Peg. A 


l Peg. B 




6642 K 


4991 K 


L/L Q 


3.041 


0.3148 


R/Rq 


1.320 


0.7499 


Rcz 1 R* 


0.929 


0.696 


Rco 1 R* 


0.100 


0.166 


T c 


17.21 


12.19 


Pc 


86.83 


79.22 


X c 


0.699 


0.682 


Y c 


0.283 


0.297 


[Li] 


+3.27 dex 


+2.69 dex 



Table 4. Calibration parameters of a iPeg. model ly- 
ing within the uncertainty boxes and of a calibrated solar 
model using the same physics. The initial values are la- 
beled by "o" . The uncertainties are empirically estimated 
(see text). 



i Peg. 



Sun 



CKA 


1 46+ - 20 


1.76 


OB 


1 -4fi+ - 30 

r.oo_ .i 




Y 


n O7S+0.007 


0.2722 


age 


0.00^°° My 


4.6 Gy 


Xo 
Z 


0.705 
0.017 


0.7081 
0.0197 



5.3. Numerics. 

The mo dels have been computed using the CESAM code 
(Morel 1997 ). The numerical schemes are fully implicit 
and their accuracy is first order for the time and third 
order for the space. Typically, each model is described 
by about 600 mass shells, it increases up to 2100 for the 
models used in seismological analysis. 



6. Results 

6.1. Age and chemical composition 

In the HR diagram Fig. |l] (bottom left) shows that the lo- 
cus of i Peg. A belongs to its uncertainty domain for evolu- 
tionary times in the interval 30 My Ss t cv A r3 500 My. Ef- 
fective temperature and luminosity are very close to their 
observed values for t cv A — 45 My. The locus of iPeg.B 
(Fig. 0, bottom right) belongs to its uncertainty domain 
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CN 


— 1 1 1 1 1 1 1 1 


/ 


_J 


/ 

i Peg A 


/ 


_J 
o 






CP 


r-it ^^-^ / l Peg/ 


B 


o 

1 






o 


i i i i i i i i 





3.8 



3.7 



3.6 



_og 



10 eff 




3.82 



3.81 



3.8 



Log 1 



3.72 



i Peg B 




3.68 



Log 1 



Fig. 1. Top: Evolutionary tracks in the HR diagram for t Peg. A and B. Bottom: enlargements around the observational 
points (full square) with the uncertainty domains. The stellar evolution sequences are initialized (stick mark *) on the 
pre-main sequence soon after the deuteron ignition. Along the stellar tracks the time intervals between the marks are 
respectively 10MY (+), 100 MY (•), 500 MY (full triangle), 1500 MY (*). 



for times 45 My t cvB 3 Gy. The distances between 
the loci of theoretical effective temperatures and luminosi- 
ties and their observed values are smaller for evolutionary 
times close to t cv b ~ 55 My and t cv b ~ 1 Gy. 

According to our definition in Sect. [|, the ZAMS oc- 
curs as soon as the nuclear dominates the gravitation as 
the primary source of energy. That occurs at evolutionary 
times t zams a — 36 My and i zams B — 70 My for i Peg. A and 
B, respectively. Therefore, the age of the system i Peg. is 
between ~ ZAMS and ~ 0.5 Gy. For the calibrated mod- 



els, within our definition of the theoretical zero age main 
sequence, the age of the binary system is about zero i.e. 
i Peg. A and i Peg. B are inhomogeneous ZAMS stars. The 
arguments in favor of this calibration are the following: 

— the low amount of lithium depletion observed in the 
secondary is in agreement with the observations, 

— in the HR diagram the loci of theoretical effective tem- 
peratures and luminosities are close to the observa- 
tions at time t GV « 56 My for primary and secondary, 
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Table 5. Theoretical global characteristics of the low de- 
gree p-mode spectrum of the star tPeg. A. All the quan- 
tities are given in /zHz (see text). 





age (My) 




Au 


So,2 


AO 


50 


2234.91 


102.0 


11.76 


Al 


150 


2222.57 


101.50 


11.47 


A2 


500 


2106.61 


96.30 


10.61 



while the differences are significantly larger as soon as 
-O.lGy << cv <0.5Gy, 
— the fact that the method of steepest descent converges 
simultaneously for primary and secondary models re- 
veals a stable solution. 

Table || shows the characteristics of l Peg. A and B cal- 
ibrated models. The primary and the secondary have a 
convective core. The lithium depletion at the surface of 
both components is close to the observed values. 

Table || shows the calibration parameters obtained for 
(.Peg. and solar models computed with the same physics. 
The uncertainties are derived empirically by changing the 
calibration parameters around their values derived of the 
calibrated model (the observed solar metallicity of l Peg. 
is assumed to be error free). Considering first, the non 
linearity of the solution around ZAMS and second, the 
large uncertainty domains, chiefly due to the inaccurate 
value of Amy, more elaborated estimations of calibration 
parameters uncertainties via, e.g. Monte Carlo procedure, 
will need heavy computations, which are not relevant of 
the scope of this paper. 

6. 2. Seismological analysis of i Peg. A 

The star l Peg. A belongs to the solar-like stars region of 
the HR diagram. The oscillations of such stars may be 
stochastically excited by the convection, as it occurs in the 
Sun. The amplitudes of solar-like oscillations have been es- 
timated by Houdek ( |1996| ) for stars of different masses and 
ages. Here, we consider the oscillation spectrum in the fre- 
quency range where the predicted amplitude is maximum, 
i.e. between 1 and 3 mHz. 

The oscillations of a non rotating star are characterized 
by two numbers: the degree t which is inversely propor- 
tional to the horizontal wavelength and the radial order 



n (Christensen-Dalsgaard & Berthomieu 1991). Here the 
rotation is small, hence neglected. The frequencies which 
may be observed, according to Houdek, are on the low- 
degree (£=0 to 3) high-frequency p-modes (i.e. large ra- 
dial order p-modes). The frequency spectrum is character- 
ized by a few global quantities, Avq which describes the 



quasi-equidistance of the frequencies in the spectrum and 
6i t i+2 the small frequency separations between frequencies 
of modes with degrees of same parity and consecutive ra- 
dial orders v n ^ and z/„-i,^+2- They can be derived from the 
following simplified polynomial expression (Berthomieu ct 
al. p993h: 



v n j = m + Av Q (n + - - n ). 

The quantity vqi varies little with the degree I. We have 
taken hq= 21, approximately in the middle of the range of 
the expected excited frequencies (i.e. 1 — 3 mHz) according 
to Houdek (loc.cit.). The large frequency separation A^o 
is the distance between peaks of consecutive radial order 
for a given degree. It corresponds to the sound travel time 
across a stellar diameter and is mainly sensitive to the 
outer layers of the star. The small frequency separations 
around n = no, defined by: 



1,1+2 



are very sensitive to the core of the star. 

Table H shows the quantities Vo,0i Avq and <5o,2 ("o = 
21) which have been computed for three models AO, Al, 
A2 of l Peg. A in the observational uncertainty HR box, at 
different evolutionary times t cv — 50, 100, 500 My. These 
characteristic quantities decrease with the age. voo and 
Avq vary proportionally to R~^, as expected, since the 
structure of the three models are very similar. The large 
separation Av$ which is the first quantity to be derived 
by observations varies significantly with the age of the 
star. The relative variation of <5o,2 is larger, due to the 
differences of structure in the core, but more difficult to 
be observed. 

7. Discussion and conclusions 

With the individual masses d erived from spectroscopic 
and visual orbits (Boden et al. 1998 ), we have attempted 
to calibrate the close binary system tPeg. We have 
adopted effective temperatures and luminosities as de- 
rived from photometric data with the K bolometric cor- 
rection and the empirical scale of temperature of Alonso 
et al. (1995, 199(i| ). We have fixed the initial metallicity of 
primary and secondary to the nowadays observed (solar) 
value. The evolutionary sequences include the pre-main 
sequence. Microscopic diffusion and gravitational settling 
of chemical species are allowed for. To fit the observed 
lithium depletion in tPeg.B we have used a mildly mix- 
ing of convection zones in the course of the pre-main se- 
quence. The convective cores are extended by 0.2 H p , the 
standard amount of overshoot. 

We have obtained adjustments within the uncer- 
tainty domains for effective temperatures and luminosi- 
ties for both components with times of evolution within 
~ 40 My ^ t cv ^ 0.5 Gy. Owing to the lack of lithium de- 
pletion observed in tPeg.B, a solution is obtained for the 
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evolutionary time t L p cg . ~ 56 My. The age of the system 
is then close to zero and we argue that the components 
of the binary system i Peg. are non-homogeneous ZAMS 
stars. 

Separate values are found for the mixing-length param- 
eter of each component. Even for t Peg. B, of mass smaller 
than solar, Table || shows significant differences between 
the convection parameter derived for each component and 
those of the calibrated solar model computed with the 
same physics. This result appears in disagreement with the 
universality of the solar convection parameter for stellar 
masses between 0.6Af Q ^ M* I.OMq with metallicity 



close to solar (Fernandes et al. 1998). Nevertheless, owing 
to the large observational uncertainty domains, our result 
does not totally rule out the paradigm of the universality 
of the solar convection parameter. 

We find a larger mixing-length parameter for i Peg. A 
(T cff = 6642 K) than for the cooler tP eg.B (T cff = 
4991 K). The calculations of Ludwig et al. (1999) seem to 
predict exactly the opposite. We do not elucidate the rea- 
sons of this difference. We are aware that our models are 
not strongly constrained by the observations. The largest 
uncertainty comes from the difficulty to estimate the mag- 
nitude difference in the V band, a fundamental quantity 
for the determination of luminosities and effective tem- 
peratures based on photometric data. Improvements of 
the calibration of iPeg. will necessitate a more accurate 
value for the V magnitude difference. Chemical compo- 
sition determinations and estimates of effective tempera- 
tures and gravities by spectroscopic detailed analysis of 
each component are perhaps realistic with modern tech- 
nology (C. van 't Veer, private communication) and are 
urgently needed. 

We must also emphasize the fact that i Peg. is proba- 
bly tidally locked or, at least, recently unlocked and then 
it maybe a chromospherically active binary. Therefore to 
model such a pair as two standard isolated single stars 
is perhaps too simplistic. In some years from now, one 
can expect that asteroseismology from space e.g. COROT 



mission (Baglin 1998) will provide more constraints on 
stellar models as it is done nowadays with the Sun. This 
preliminary work also shows the powerful interest of in- 
terferometric and spectroscopic binaries such as (-Peg. as 
tests of transport process theories in stellar modeling. 
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